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Abstract 

We present a set of formulae to extract the gluon distribution function from 
the deep inelastic structure function F2 and its derivative dF2/dlnQ^ at small 
X in the leading and next-to-leading order of perturbation theory. The detailed 
analysis is given for new HERA data. The values of the gluon distribution are 
found in the range 10'^ <x< lO^^ at Q"^ = 20 GeV^ 
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The knowledge of the DIS structure functions at small values of the Bjorken scaling 
variable x is interesting for understanding the inner structure of hadrons. Of great rele- 
vance is the determination of the gluon density at low x, where gluons are expected to be 
dominant, because it could be a test of perturbative QCD or a probe of new effects, and 
also because it is the basic ingredient in many other calculations of different high energy 
hadronic processes. 

Recently two experiments working in the electron-proton collider HERA at DESY (HI 
and ZEUS) have published new data on the structure function F2 ||T|], 0. Up to now all 
the analysis performed of these data 0-0 found that the gluon distribution rises steeply 
towards low x (in the moderate range of the measurements). This behaviour has been 
recently connected within the DGLAP evolution equations with the less singular are 
found at lower values of by NMC and E665 experiments. 

We introduce the standard parameterizations of singlet quark s{x,Q'^) and gluon 
gi^XjQ"^) parton distribution functions (PDF)| (see, for example, 0) 

s{x,Q^) = AsX'\l-xY'{l + es^/x + -fsx) =x'^s{x,Q'^) 

gix,Q^) = A,x-\l-xY^{l + eg^/E + ^gx)=x-'gix,Q^), (1) 

with dependent parameters in the r.h.s.. 

Note that the behaviour of Eq. (|I|) with a Q^-independent value for 6 {6q = 6g) obeys 
the DGLAP equation when x-^ > 1 (see, for example, @ - 0). If S{Ql) = m some 
point Qo ~ IGeV'^ (see |Tl|, |]12|, 0) , then the behaviour p{x,Q'^) ~ Const {p = {s,g)) 



is not compatible with DGLAP equation and a more singular behaviour is generated. If 
we restrict the analysis to a Regge-like form of structure functions, one obtains (see P) 



p{x, Q"^) ~ X 



with next-to-leading order (NLO) 5q(Q^) 7^ ^giQ"^) intercept trajectories. 
Without any restriction the double-logarithmical behaviour, i.e. 



p{x,Q')r^exp{yS,{Q')ln^) (2) 
is generated. At NLO and for / = 3 active quarks, 5p(Q^) have the form (see 0): 



■'9 

while for f = 4: one has: 



3 81 



where t = /n(a((5o)/a(Q^)) and Z = a{Ql) — a{Q'^). 

Because we would like to get formulae to extract the gluon distribution from exper- 
imental data without theoretical restrictions, we will consider both, the Regge-like be- 
haviour of Eq. (P if x~^ 3> 1, and the non- Regge-like behaviour of Eq. (0) if S{Ql) = 0. 

^We use PDF multiplied by x (s = xS and g = xG) and neglect the nonsinglet quark distribution at 
small X 
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In this work we present a simple relation between the gluon and F2 which permits the 
extraction of the gluon directly from data with a NLO precision in perturbative QCD. 
This kind of formulae (see also the NLO result of ref. [|l^] which depends on an unknown 
phenomenological function) are useful as an alternative to the most complex analysis in- 
volved in global QCD fits. At low x the coupled integro-differential equations can be 
converted in more simple linear relations between parton densities and structure func- 
tionsQ. The method to arrive to the solution presented below is based in the replacement 
of the Mellin convolution by ordinary products developed in ref. that was already 



applied in the derivation of the leading order (LO) formula in ref. and the NLO result 
1/2 in ref. 



at S 



1. Assuming the Regge-like behaviour ior the gluon distribution and F2{x,Q'^): 
g{x, = x-'~g{x, Q'), F^^x, Q^) = x-'f{x, Q'), 
we obtain the following equation for the derivative of the SF F-Q. 



dF2{x,Q' 
dlnQ^ 



— X 
2 



(3) 



where r^p(a) are the combinations of the anomalous dimensions (AD) of Wilson operators 
7^^p = a7g)'^ + 027(1)'^ + 0{a^) and Wilson coefficients^ aEf''' + 0{a^) of the r] "moment" 
(i.e., the corresponding variables expanded from integer values of argument to non- integer 
ones) : 



(4) 



gg 



7, 



(0),r)- 



and 



p'{0, Q2) ^ -^p{x, Q^) at X = 
ax 



where e = J2i is the sum of squares of quark charges. 

For the gluon part from r.h.s of Eq.(^ with accuracy of O(x^), we have the form: 



rlt'~g{x/^sa,Q') with 



sg 



I s 



(5) 



In the quark part a similar simple form is absent because the corresponding LO anoma- 
lous dimensions, 75°-*'^^^ and 7^°''''', have opposite signs. However, within accuracy O(x^), 
it may be represented as a sum of two terms like Eq.(^), with a shift of some coefficients 
and arguments. Choosing the shifts as 1 and A we have the following representation for 
the quark part: 

cis{x,Q^) + C2s{Ax,Q^) + 0{Ax^), 

^There are analogous formulae connecting with F2 and its derivative and also for the extraction 
of gluons from (see the at LO refs.[|l^ and and at NLO refs. |Q| and respectively) 
^Hereafter contrary to the standard case we use a{Q^) = as(Q^)/47r . 
^Because we consider here F2{x,Q^) but not the singlet quark distribution 
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where 



C2 



and ci = 



A-1 ^ A-1 

We thus find the following expression from Eqs. (|5)-(|§) 

^^^^ = -I [rlfiW'aix/^s,, Q') + c,F,ix, Q') + c,F,{Ax, Q') 
because q ~ 0{a). 

From Eq. (|^) with the accuracy of Oi^x"^^^ , aAx"^^^), we obtain for gluon PDF: 



(6) 



(7) 



g{x,Q^) 



sg) 



2 ■ 



dF2{xisg.Q' 



dlnQ^ 



+ Oix'^-^^aAx'^ 



Because the value of A is arbitrary it is convenient to neglect the contribution from 
F2{x^sg/ A,Q'^). Putting formally A = oo[], one arrives to the final general formula to 
extract g{x,Q'^) in the NLO approximation: 



gix,Q^' 



sg) 



(9) 



Restricting the analysis to 0{x^ ^^ax^ '^), one can replace isg ^ i = ifg'^^^ llfg'^ 
into Eq. (1): 



gix,Q^ 



sg 



2 ■ "^^f^Q?'^ + rlt' F,{x^,Q^) + 0{x'-\ax'-^) 



(10) 



This replacement is very useful. The NLO AD 7^^'*'" are singular^ in both points, 
n = 1 and n = 0, and their presence into the arguments of Q"^) makes the numerical 
agreement between this approximate formula and the exact calculation worse (we have 
checked this point using some MRS sets of parton distributions). 

Using NLO approximation of r^^^ we easily obtain the final results for g[x, Q^): 

e \dF^{xi,Q^) 



g{x,Q 



2/ 



ae 



(0),l+(5 



Isg 



+ Tsg OL 



dlnQ^ 



a 

+ 2'1 7(0).i+^F2(a:e,Q^ 



7(°)'^+^ + 7il)'^-^^«)F2(xe,g2) + o(. 

+ 0{a^,x^-\ax^-^) 



g{x,Q') 



2/ 



1 



ae 



(0),1+S _(l),l+<5 



Isg 



+ Tsg" a 



dF2{x,Q^ 
dlnQ^ 



+ 



a { (7(o).i+5 + 7(i)'i+^a)F2(x,g2) + 0{x^- 



0{a ,x 



2 ^l-S\ 



:i2) 



^Really A ~ (x^g)-^ 

®In the case of replacement Mellin convolution by ordinary product these singularities transform to 
logarithmically increasing terms (see and [^ ) 
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where 



I sg 



{0),»? 

99 



7, 



/ ss 



Any equation from above formulae (10) may be used, because there is a strong cance- 
lation between the shifts in the arguments of the function F2 and its derivative and the 
shifts in the coefficients in front of them. 

For concrete values of 5 = 0.5 (see also [1^) and 5 = 0.3 we obtain (for f=4 and MS 
scheme): 



if 5 = 0.5 

g{x,Q')- 



105 



4eV23 X 77 a (1 + 26.93a) 



dF2mx,Q^] 
dlnQ^ 



16 /107 



g{x,Q' 



105 1 



60 



2ln2 



{l + 37.76a}F2(ffx,g^ 



1 



92e a (1 + 26.93a) 
2/n2 



16 /107 

H a 

3 ^ 60 



dF2ix,Q^) 
dlnQ'^ 

\ {l + 37.76a}F2(x,g2 



0{x^ ",ax 



2-5 
2 



+ 0{a^,x 



+ 0{x^-^) 
+ 0{a^,x^-^) 



(13) 



(14) 



if 6 = 0.3 



g{x,Q' 



0.60 



ae (1 + 52.52a) 



dF2ix/5.27,Q^ 
dlnQ^ 



+ 189a I {l-50.67a}F2(x/5.27,g2 
1 F2(x/5.27,Q2) 



+ 0{x 
+ 0{a'^,x^-\ax^-^) 



(15) 



g{x,Q' 



0.98 



ae (1 + 52.52a) 



dF2{x,Q'^ 
dlnQ^ 



r {l- 50.67a}F2(x,g2) + 0{x' , 
■ 1 F2(x,Q2) + 0(a2,xi-^) 



(16) 



2. Assuming the non-Regge-like behaviour foT the gluon distribution and F2{x,Q 



2\ 



exp(i.MJQ2);^l) exp (i A(Q2)/rii) ^ 

(27r53(Q2)/ni)'/' (27r5,(g2)/ni)^/' 
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we obtain the following equation for the derivative of the SF F^: 

c/F2(x,g2) 



1 ^ exp (kxIbAQ^^ln-) , , ^ . 

^ v=s,9 {2Trdp{Q^)ln^) 



(17) 



where rl (a) can be obtained from corresponding functions rl^^{a) replacing the singular 



sp 



sp 



term 1/5 at 5 ^ by another term 1/6: 

1 



1 5^0 1 

6 ^ 6 



ln{l/x) 



E 

m=l 



1 X 3 X ... X (2m- 1) 



The singular term appears only in the NLO part of the AD 7^p'*'^''"'^ in Eq. 
replacement (|18D corresponds to the following transformation: 



(l),l+<5 ^ ^(1),1 
I sp I sp 



^ / sp I sp 



^W.il , ;^(i),i 





The 



(19) 



where 7^^'*'^ and 7(^)'^+'^ are the coefficients corresponding to singular and regular parts of 



I sp 



7sp'''^^^5 respectively. 

We restrict here our calculations to 0{x) because at O(x^) one obtains an additional 
f cLct or ' 

exp {l^6,{Q^)ln^) {27r6,{Qynlf' 

in front of the function F2 and its derivative. This factor complicates very much the final 
formulae. 

Repeating the analysis of the previous section step by step using the replacement (p!9D, 
we get (for f=4): 



g{x,Q^ 



3 1 



1 



4ea(l + 26a[l/5- §]) 



dF2{x,Q^) 
dlnQ^ 



a 



'{203 - 61/5}F2(x,Q2 



O(x^) 
0{a^x'] 



(20) 



because 7^°''^ = 0. 



3. In the case x^^ ^ Const our formula (|T2D at accuracy 0{x^~^) coincides with 
the corresponding Ellis-Kunszt-Levin result from ref. |T0[. By other part, for 6 = 1 one 
arrives to 



gix,Q^ 



1 



Aea (1 + 



32 
+ — a 
9 



dF2{x/2,Q^) 
dlnQ"^ 

26231 1181/]«}F2(x/2,Q2) + 0( 



1728 



576 



x'^'^ax^-^) 



Oi^a^.x^ ^,ax 



l-5\ 



(21) 



^Using a lower approximation 0{x) is not very exact, because in this case F2 and the gluon distribution 
can contain an additional factor in the form of a serie 1 + '^/^{^ / Sp / ln{l / x))'^ , which is determined by 
boundary conditions (see discussion in Ref.[|l2j). We will not consider the appearance of this factor in 
our analysis 
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that coincides with Prytz results |TH] in the LO approximation, when we neglect the 
contributions ~ F2{x,Q'^). Both formulae, ((|2T|) and one from |]16|), are similar to ours 
in the NLO case, too. Certainly, the value 6 = 1 lies outside the more standard predicted 
range < 6 < 1/2, however in the case of large 6 {6 > 0.25) the final formula (|10|) depends 
very slowly from the concrete value of 6. This is due to the strong cancelation between 
the shifts in the arguments and in the coefficients in front of the functions. 

Equations (IT^) and (12^) can be combined in a more general formula valid for any value 
of 5: 



g{x,Q'^ 



2/ 



1 



dF2{x,Q^) 



+ 



a 



(0),l+<5 



dlnQ^ 



+ 0{x^-^) 
+ 0{a^,x^-^) 



(22) 



where J^^^'^'^^ coincides with 7^1,^'"'^'*''^ with the replacement: 



1 

6 



dy giy,Q^ 



y 9{x,Q^ 



(23) 



In the cases x 3> Const and 5^0 the r.h.s. of (p3D leads to 1/5 and 1/S, respec- 
tively. 



4. In Fig. 1 it is shown the accuracy of Eqs. (|T3|), (|TBp (both 0(a;^,x^ ^,ax^ ^)) 
and ( pOD {0{a'^,x^)) in the reconstruction of various gluon distributions from MRS sets 
at Q2=20 GeV^. We have chosen for this test MRS(Do) (5=0), MRS(D_) (5=0.5) and 
MRS(G) (5=0.3) as three representative densities (see ref. |^ and references therein). It 
can be observed in Fig. la that using the formula with 5 = 0.5 one gets the best agreement 
with the input parameterization (less than 1 %) in the case of MRS(D_) set; for MRS(G) 
the reconstruction is still good (less than 10 %), but for MRS(Do) the deviation reach a 
30 % at low X. 

In Fig. lb the degree of accuracy of the reconstruction formula with 5 = 0.3 can be 
observed. Here one should expect the set MRS (G) to give also a very good (~ 1% level) 
agreement, however this is not the case because set (G) distinguishes the exponents of the 
sea-quark part 5^ ~ from the gluon density (5 = 0.3). Thus, Eq. (12) might be slightly 
modified to treat this case. 

Figs. Ic and Id deal with the case 5 = 0. As in Fig. la, one can observe a very good 
accuracy in the reconstruction when coincides with that of the test parameterization 
(4 GeV"^ for MRS set). Notice also the lost of accuracy at high x due to the importance 
of the 0{x) terms neglected in Eq. (16) 

With the help of Eq. (|1^) we have extracted the gluon distribution from HERA data, 
using the slopes dF2/dlnQ^ determined in ref. and ref. [Q. When HI data are used 
the value of F2 in Eq. (|14|) was directly taken from the parameterization given by HI in 
ref. With ZEUS data we substitute directly the F2 values presented in table 1 of their 
ref. IQ. We have checked that the use of the HI parameterization for F2 when dealing 
with ZEUS data, does not change significantly the xG{x, Q"^) result. 
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Figures 2a and 2b shows the extracted values of the gluon distribution. It can be ob- 
served that the agreement within the errors between the bands, generated from a global 
fit to data, the parameterization MRS(G), and the extracted points is excellent. 

5. In conclusion, a set of new formulae connecting the gluon density with F2 at low x 
have been presented. They work fairly well for singular type gluons {6 ~ 0.5 — 0.3) and 
for the non-singular case (5 ~ 0). We have reproduced previous results of Prytz [0 using 
6=1, and of Ellis-Kunszt-Levin |]T^ with accuracy 0{x^~^). 

We have found that for singular type of gluons the results do not depend practically 
on the concrete value of the slope S: there is a cancelation between the changes in the 
arguments and in coefficients in front of the functions. However, when 5 — the co- 
efficients in front of dF2{x , Q"^) / dlnQ"^ and F2{x,Q^) have singularities leading to terms 
~ ln{l / x)^ Consequently, before to apply these formulae, some fit (maybe quite crude) 
of experimental data is necessary to verify the type of F2{x) asymptotic at x — * 0. 

The formulae were used to generate the gluon distribution that agree with the rise 
observed by HI and ZEUS experiments. Further work is in progress in order to obtain 
similar expressions connecting F^, F2 and the derivative of F2. 
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Figure captions 

Figure 1: Relative difference between the reconstructed gluon distribution using formulae 
in text and different input parameterizations 

Figure 2: The gluon density. The points were extracted from Eq. (^^ using HI (Fig. 2a) 
and ZEUS (Fig. 2b) data. The dashed curves shows the limits of the error band taken 
from Fig. 3a of paper [Q and Fig. 4 in ref. [Q which represents the uncertainty from a 
NLO fit. Solid line is the gluon density from set MRS(G) g] 
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